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Abstract 

In this paper we discuss a new way to derive neutrino mixing patterns, which originates from 
the idea proposed in a recent article by Hernandez and Smirnov. Its applications to various cases 
are discussed. We first present the complete set of possible mixing patterns for the minimal case 
where unbroken residual symmetries of the Major ana neutrino and left-handed charged- lepton 
mass matrices obey some general assumptions that are also satisfied by many models based on 
discrete symmetries. We find that they are either well-known mixing patterns or phenomenolog- 
ically disfavored ones. It shows clearly that, for full-mixing matrices to fit the mixing data with 
small or negligible corrections, it is necessary go beyond the minimal scenario. We present an 
explicit formalism for a rather general non-minimal case. Some applications and phenomenological 
implications are discussed. Several new mixing patterns are derived. 
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I. INTRODUCTION 



The discovery of neutrino oscillations provides a clear evidence for physics beyond the 



Standard Model (see, e.g. [l, 2] and for recent global fits, see 
neutrino mixing determined by underlying discrete flavor symmetries 
on this approach have been discussed extensively (see, e.g. UQ for some recent works) 
and they often lead to intriguing mixing patterns, including the well-known tri-bimaximal 
mixing (TBM) 32] and bimaximal mixing (BM) 33| . However, a challenge is posed by recent 



-|5j). In this work we study 



aa. 



Models based 



34-38 



because many models discussed 



experimental results on the reactor mixing angle 
previously did not predict the measured value of #i 3 (see, e.g. |39j or [40j for a survey of 
model predictions). Although good agreement with the data can be achieved by introducing 
corrections or other means (see, e.g. 4l|-|44l|). it is still worthwhile to seek for symmetries 
and mixing patterns that are in good agreement with the data without introducing large 
corrections which may badly break the (would-be) residual symmetries, e.g. the THF mixing 



45 



47]. 



or bi-trimaximal mixing discussed in recent articles 

In this paper we focus on Majorana neutrinos. The mass matrix of Majorana neutrinos 
has a Z 2 x Z 2 symmetry and the left-handed charged lepton mass matrix has a U(l) x U(l) 



symmetry 



48 



49l | . if they are required to belong to SU(3). When restricted to discrete 



symmetries, the later is usually reduced to a Z m or other finite symmetry belonging to the 
U(l) x U(l) (for some examples, see, e.g. 45]). In models based on discrete symmetries, 
these symmetries may coincide with the symmetries preserved by relevant mass terms after 
the breaking of flavor symmetry 1 , which is the case that will be considered in this paper. 



We follow 



50] and carry out a model independent study of neutrino mixing. The basic 



assumption we adopt is that the residual symmetries of the mass matrices belong to a 
underlying discrete flavor symmetry. Our main concern is the general results that can be 
drawn from this assumption. 

We find that the idea of 50] leads to a new way to derive mixing matrices, which is 
simpler than the group theoretical method. It does not require a flavor group to be specified 
in advance. On the contrary, it can provide the necessary conditions form which one can find 
the flavor symmetry corresponding to a mixing matrix or mixing-matrix column derived by 



1 Models with completely broken symmetries can be constructed (e.g. [ssj]), which are beyond the scope of 
this work and hence will not be considered here. 
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this method. If the flavor symmetry is known, it can serve as a quick way to derive mixing 
matrices or elements and cross check with group theoretical calculations. From the equation 
constraining the residual symmetries, an interesting result we find is that, under some general 



and plausible assumptions which are satisfied by many models and adopted in |50j, one can 
prove rigorously that phenomenologically viable mixing patterns that discrete symmetries 
can lead to include only TBM, BM and the golden ratio mixings. It also provides a different 
way to approach some known results. For example, one can show that the minimal horizontal 
symmetry that can produce TBM as a full-mixing matrix is S4, which was first pointed out 



in 



5l|- 



As suggested by the no-go result mentioned above, those assumptions may have to be 
relaxed, and hence new mixing patterns can be obtained. Explicit formalism for a rather 
general non-minimal case are also provided in this paper. Several interesting applications 
are discussed and new mixing matrices are presented. Their phenomenological implications 
are also discussed briefly. We find that, although not impossible at all, it is still difficult to 
find a mixing matrix that can produce all the mixing data. This result agrees with [53 1 and 
54j in which the results of the searches for groups up to certain large orders are reported. 

Although our approach is based on that o f |50l| , the concerns and discussions differ quite 



a lot. We also note that, in a recent work [52j, authors of [50] also generalize the results 
of their previous work. But the subjects and the formalism are still quite different. More 



attention is paid to partial-mixings and phenomenology in [52j and their discussions are more 
group theoretical oriented. Most of our attention is concentrated on full-mixing patterns. A 
general result about the minimal scenario of [50] is presented and discussed in this paper. 
Besides the formalism that can be used in general cases, we also provide a simple formalism 
for a rather general non-minimal scenario from which most known mixing patterns, including 
those discussed in {52] or 45], can be derived. It also provides a different way to understand 
the permutation feature of some mixing patterns. Several new mixing patterns are also 
derived and presented. 

The paper is organized as follows. In Sec. 2, we review briefly the idea proposed in 



50j. For later convenience, the formalism for squared mixing elements, which can be used 



in general cases, is presented in this section. In Sec. 3, the attention is focused on the 
minimal scenario. A general result concerning possible mixing patterns is presented and 
discussed. Then, in Sec. 4, the consequences of relaxing some of the assumptions adopted 
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in the minimal scenario and the relevant phenomenology are discussed. We summarize in 
Sec. 5. 



II. NEUTRINO MIXING AND DISCRETE SYMMETRIES 

We first set up our notations. We denote the generator of the symmetry of the left-handed 
charged-lepton mass matrix (i.e. M\M x where M { is the mass matrix of charged leptons) 
by T and those of the Majorana neutrino mass matrix by Si where i — 1, 2 or 3. The 
flavor symmetry is denoted by Gf and the symmetries generated by T and Si by G e and 
G u , respectively. In this work, we always assume that G v = Z 2 X Z 2 and T generates a Z m 
symmetry. We also assume that T and Si belong to SU(3) and, for simplicity, Tr[T] is real 
in most of our discussions. Some comments on these assumptions are given at the end of 
this section. 

In the basis where the charged lepton mass matrix is diagonal, T can be written as a 
diagonal matrix denoted by T^. In the case where Tr[T] is real, Td can be written as one of 
the following three matrices, 

T e = T 1 = diag{l, e 27Tik/m , e ~ 2 ™ fc / m }, 

T M = T 2 = diag{e 2 " fc/m , 1, e " 2 ™ fc / m }, (1) 
T T = T 3 = diag{e 27rik/m , e~ 2nik/m , 1}, 

since it belongs to SU(3). Therefore, 

Tr[T] = 1 + 2cos2vr£;/m. (2) 

when Tr[T] is real. In addition, we assume that T a (a = e, fi, or r) is non-degenerate or 
m > 3, otherwise charged lepton masses cannot be guaranteed to be non-degenerate and 



ambiguities in mixings may arise 491 ] . Further discussions can be found at the end of this 
section. 

The key assumption on T and Si is that they belong to Gf, a finite flavor symmetry. It 



then follows that (SiT) Pi = I for some Pi which leads to 50| 



{SiT) Pi = (UpMNsStul MNS T a y> = I (3) 



2 Other relations are possible if T or Si appears in Gf indirectly. 
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where Upmns is the lepton mixing matrix and Sf are diagonal matrices given by 

S? = diag{l,-1,-1}, S* = diag{-l,l,-l}, Sl = SfSl (4) 
Then denote (S l T)- 1 = T^Si by W % . One finds that 

{W i r=[{S,{T)- 1 ] Pl =l. (5) 
Our basic assumption on the residual symmetries can then be summarized by the relations 3 

S 2 = T m = fyy^pi = ^rp^ = j ^ 



The characteristic equation of Wi can be written as [50 1 



A 3 + WiX 2 - w*X -1 = (7) 

where Wi = — Tr[Wi]. Because of Eq.([5]), one has X Pi = 1 and from the equation above it 
follows that if Wi is real, the eigenvalues of Wi are given by A = 1, e 2nmi ^ Pi , or e - 27n? Wp* an d 
hence Wi = —1 — 2coa2i:n i /p i . 

Until now we have followed 1 501 . Below we shall discuss first the general situation. Later 

n 

in the next section we return back to the minimal scenario discussed in [50] . Since 

Tr[W,] = TriiSiT)' 1 } = Tr[(^T)t] = (Tr[5,T])* , 

one finds that 

w* = - (Ti[W t ]Y = -Tr[UpMNsStU PMNS T d ] = Tr[T d ] - 2Tr [AiT d ] (8) 
where Aj are defined as 



Ai — diag{| (UpMNs)u | 2 ' I (UpMNs)2i | 2 > I {UpMNs) Si | 2 }- 

Eq. ([8]) is the most general formalism for mixing elements in this framework. As will be shown 
later, it can be written in more convenient and explicit forms without matrix operations. 



3 It is pointed out in 50 1 that these relations define the von Dyck group D(2,m,p) and the only finite von 
Dyck groups are S 3 , S4, A4, A 5 and D(2,2,p), the dihedral groups. However, other finite groups can be 
obtained by imposing additional relations. Therefore, our formalism can be applied to those groups as 
well. More relevant discussions can be found at the end of Sec. 3. 
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From Eq.®, an important condition can be derived. Summing over % leads to 

^< = -^(Tr[^]r = Tr[T d ] (9) 



i=i i=i 

which follows from Yli=i 8f = ~^ an d the unitarity of Upmns- This equation will be referred 
as unitarity condition in this paper. Since this condition always holds, it must be obeyed by 
any combination of Si and T that generates a full- mixing matrix. As Eq.flSJ), it can also be 
written in more explicit forms. Eq.® and (Q are the starting point of our discussion that 
follows. 

When Tr[T] is real, is given by one of the T matrices in Eq.([T]). Then from Eq.® it 
follows that 

ReH = 1 - 2| (U PMNS ) ai | 2 (1 - cos2nk/m) (10) 
lm[wi] = 2 | (Upmns)^ | 2 - I (UpMNs) 7 i I 2 sm2irk/m (11) 

where /?, 7 7^ a and /? < 7. Solutions to Eq. (|T0|) and (|TT|) should respect the conditions 

0<|(f/pM7V5) pi | 2 <l, p = a, /3, or 7 (12) 

which will be referred as reality conditions. Explicit expressions for | (Upmns) ~ | 2 wm be 
given later. Note that the value of a depends on which T matrix given in Eq. ([1]) is used in 
Eq.([81). Using a different T matrix but keeping m and k fixed results in a reordering of the 
elements in the mixing vector by which we mean a column of the mixing matrix. One may 
use any T matrix given in Eq.flTJ as long as the results are consistent with the experimental 
data. 

Together with the unitarity and the reality conditions, the above formalism provide a 



simple way to derive mixing matrix or mixing elements. For example, it was shown in 45] 
that the group PSL(2, Z-j) can lead to a mixing matrix given by 

1 V2 1 



I^PAfiVsll — ^ 



^(3-v/T) 1 ^(3 + y/f) J 



(13) 



4 Following 45[, we use ||J7pMiVs|| to denote the matrix with every entry being the absolute value of the 
corresponding one in the mixing matrix Upmns- 
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The representation matrices of PSL(2, Zj) are given by 



2 

7f 



^ s l s 2 S3 ^ 



52 -S3 Sl 

53 Si S 2 



( p 4iri/7 



G 



2 

7f 



V 










e 2ni/7 








e 8ni/7 



(14) 



where sj~ = sin/c7r/7. Si, S2 and T can be chosen as F,G 2 FG 3 FG and G 3 F, respectively. 
The ||?7pAfJvs|| matrix given in Eq.f JT3|) can be obtained by diagonalizing Si and T. The 
formalism given above provides a different but more efficient way to derive it. More explicit 
formalism and examples will be given later. 

Now we have all the necessary ingredients, but before proceeding, there are a few com- 
ments we would like to make: 

1. Since T G SU(3), if Tr[T] is not real, then it does not have a +1 eigenvalue and hence 
no vacuum alignment can break Gf into the Z m symmetry generated by T. It can 



be arranged that residual symmetries are preserved indirectly |46|, l55[ , but it sounds 
more natural if they emerge directly and thus, except in an example given in Sec. 4, 
we assume in most of our discussions that Tr[T] is real. 

2. When degeneracy occurs among the eigenvalues of T, whether Tr[T] is real or not, 
lepton mixing matrix can not be determined unambiguously and the solutions to Eq.flH]) 
are not unique. Nevertheless, the vanishing of Im(Tr[T]) assures that for any m > 2 
the eigenvalues of T are non-degenerate and the mixing matrix can be determined 
unambiguously. As to the case where m — 2, to eliminate the ambiguity, one may 
enlarge the Z 2 symmetry generated by T to a larger one, e.g. a Z 2 x Z 2 . An example 
is given in footnote 5. Since we are interested in general cases, we will not consider 
this particular case further. 

3. Unlike the assumption on Tr[T], Tr[Wi] are assumed to be real mostly for calculational 
simplicity. In many cases, including the PSL(2, Z 7 ) example discussed above, Tr[Wj] 
are not real. But we should mention that, this assumption is indeed satisfied by many 
models leading to well-known mixing patterns (e.g. TBM, BM and the golden ratio 
mixings). In addition, it can be treated as a reasonable phenomenological assumption. 
The reason is that, the equalities between the absolute values of mixing elements 
are still phenomenologically viable and from Eq.f llip it follows that when Eq. (ll6j) is 



satisfied, Im[u>i] vanishes automatically. Detailed discussion about the case where 
Tr[Wj] are not real will be presented in Sec. 4. 

4. From experimental data it follows that in the minimal scenario sin #13 can only result 
from Eq. (|15p (see below) and hence, to accommodate measured sin# 13 ~ 0.15, 1 — Wi 
must be almost vanishing because sin 2 9\s ~ 0.02, which requires a large (recall that 
Wi = — 1 — 2cos27rnj/pi) such that ni/pi can be close to 1/2. Numerical computation 
also shows that to fit the data, Pi should be larger than 10. Large Pi corresponds to a 
S{T of large order. Therefore, the larger the Pi is, the less the chance for small groups 
to accommodate the data. Allowing corrections to sin #13 may significantly lower the 
requirement, but introducing large corrections may also require the price of weakening 
the role of symmetry to be paid. More discussion can be found in Sec. 4. 

5. If all Si G Gf and m > 3, then the matrix H^pmtvsII is completely determined by 
Si and T and is referred as full- mixing matrix. In the case where both Tr[T] and 
Tr[S"jT] are real, one can show that Eq. (fT6|) (see below) always lead to a vanishing 
mixing element or a maximum Dirac phase in an appropriate parameterization (56 ]. 
The vanishing element implies that Wi = 1 and thus from Eq. (J7j) , one finds that Pi = 2 
and (SiT) 2 = 1, i.e. one S^T must be of order 2, which can also be seen from the result 
presented in the next section. Note that Pi cannot completely determined in this way. 



III. MIXING PATTERNS IN THE MINIMAL SCENARIO 



In this section, we focus on the minimal scenario where both Tr[Wi] and Tr[T] are real. 



In this case, the solution to Eq.f JTOj) and f|TT]) can be written as [50|: 



2 l-Wj 1 + TrtWj 

I [UPMNSjai I - „ • 2 kn = „ • 2 kn t ib > 

4 sin — 4 sin — 

m m 

I {UpMNs)pi I 2 = I {UpMNs)-,i I 2 = 2 ( X ~ I ( U PMNs) a i H • (16) 

We also assume that all ^ G Gf. Hence pi and m are strongly constrained by the 
unitarity condition given by Eq.flUJ). As discussed below Eq.flT]), if TrfWj is real, then 

ry ■ 

Tr[Wi] = 1 + 2 cos— 2tt (17) 

Pi 
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where pi is the order of W{. Note that the greatest common divisor of rii and pt, gcd(nj, pi 
1. Now from Eq.([2]) and (Q, one has 



, 3 3 

k \ _ , \ m 



1 + 2 cos — 2tt = -^2 Tr t^] = -3 - 2 cos ^ 27r ( 18 ) 
which can be written as 



m * — ' * — ' va 

i=l i=l ^ 



4 

n 



^2cos-^2vr = -4 (19) 

3=1 Pj 

where is replaced by j?4 and m by n.4. Without loss of generality, we require that < 
n j/Pj < 1/2- Eq. (fl9|) is a necessary condition that pj, the orders of T and VFj (or S^T) 
must obey. Together with rii, they can be used to construct explicitly the (minimal) group 
corresponding to a solution. Its relation with the group theoretical method will be discussed 
shortly. 

We find that a c omp lete and rigorous solution to Eq. ffl9|) can be derived by using algebraic 



number theory 58|,|59j], which can translate Eq. (ll9j) into a much simpler arithmetic equation. 



Although a rigorous derivation is interesting by its own right, it is somewhat lengthy because 
some mathematical concepts and results need to be introduced. In addition, the solutions 
given below can be verified numerically by evaluating Eq. ( fl9l) continuously until pt reach a 
desired large value. Therefore, in order to concentrate on physical discussion, we will present 
the detailed mathematical derivation elsewhere. 

Since can be easily found from Eq. (ll9j) when pi are given, below we denote the solution 
to Eq. (fT9|) by {pj}- Although expected, it is still a little bit surprising to find that, besides 
obviously disfavored solutions, i.e. {1,2,2,2} and {2,2,p 3 ,p 4 }, Eq.(fl9l) admits only three 
other solutions: {3,3,3,3}, {2,3,3,4} and {2,3,5,5}. Note that ps and p^ in the second 
solution satisfy n 3 /p 3 + n^/p^ = 1/2 where n 3 and n 4 are arbitrary integers satisfying 
< rij/pj < 1/2 and gcd(rij,pj) = 1. 

To find the ||?7pmjvs'|| matrix corresponding to a solution, one first assign one pj to m 
and the others to ^1,2,3 i n Eq. ([T7|) . It is obvious that the first solution {1,2,2,2} does not 
respect the assumption that m > 3 and hence it cannot lead to an unambiguous mixing 
matrix. The second solution {2,2,^3,^4} is also phenomenologically disfavored because of 
the same reason, which requires that any pj = 2 can only be assigned to the order of a S^T. 
From Eq. (fl5|) and ( ITT)) it follows that U ai = if Pi = 2. However, p 1>2 = 2 result in two 
vanishing elements which lead to two vanishing mixing angles and hence this solution is not 
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phenomenologically favorable. 

The third solution, {pi, m} = {3, 3, 3, 3} leads to an interesting ||?7pmjvs|| matrix in which 
all the elements equal to 1/Vo, but it is also not phenomenologically favorable because of 
its large deviation from the experimental data. Note that it leads to a maximum Dirac CP 
phase. 

Therefore, we are left with the last two solutions. From Eqs.( ll5fH7|) . it is easy to show 
that the assignment {pi,p 2 ,P3,w} = {3,3,2,4}, {3,4,2,3}, {3,5,2,5}, or {5,5,2,3} leads 
to TBM, BM, or the golden ratio mixings, respectively 5 . As an example, let {pi,P2,P3, m } — 
{3, 5, 2, 5} and {ni, n 2 , n 3 , k} = {1, 2, 1, 1}. From Eq.fJTTJ) and Tr[Wi] = —Wi, one has 

Wl = 0, W2 = -l-2 cos[4tt/5] = (V5- l)/2, w 3 = 1. 

Then from Eq.flTS]) and 4 sin 2 tt/5 = 10/(5 + it follows that 

| (UpMNs) a l I 2 = 1Q - ' I (UpMNs) a2 ? = 1 ' (U PM Ns) a 3 P = 



Other elements can be obtained from 
matrix (for explicit expression, see 



45 



BP and the resulting Upmns is a golden ratio mixing 



60j). Note that, one must set a = e = 1 for phe- 
nomenological reasons. Also note that, mixing matrices obtained by exchanging pi or using 
a different value for a are just the original mixing matrix with it rows and columns being re- 
ordered. But for the assignments discussed above reordering is obviously phenomenologically 
unacceptable. 

This example shows clearly the difference between the method discussed here and the 
group theoretical method. The later requires the knowledge of the flavor group and its 
representation. Here what we need is a set of Pi and that satisfy the unitarity condition and 
no prior knowledge of the flavor group is needed, which can be determined by the solution. 
For instance, the example discussed above corresponds to the solution {p^ = {2,3,5,5} 
which indicates that the flavor group must contain Z 2 , Z 3 , Z 5 , and Z 2 x Z 2 subgroups 
and its order must be a multiple of 60. It is not hard to find that the minimal finite group 



5 As discussed in the last section, we require that m ^ 2. If T 2 = I, to fix the mixings, one may enlarge 
the Z2 symmetry generated by T to a larger one. For example, if {3,5,5,2} is assigned to {pi,m} and 
G e is chosen to be Z2 x Z2, then G e contains two Z2 generators, each of which determines a row of U. 
Therefore, U can also be completely determined in this case. The corresponding mixing matrix can be 



found in 



45|. 
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satisfying this condition is A5. One can use pi and n« given above to verify that this solution 
can indeed be realized by A5. Similarly, from the solution leading to TBM one can show 
that S4 is the minimal horizontal symmetry that can produce TBM as a full-mixing matrix, 



which was first pointed out in 



51] . In addition, one may also use Pi and rii to construct the 



group. In summary, the solution associated with a mixing matrix can be regarded as the 
necessary conditions or minimal requirements for a flavor symmetry to produce the mixing 
matrix. Although necessary conditions are not sufficient to establish the existence of a finite 
group, but so far no counter-example has been found. 

In the discussion above it is assumed that all Si and T belong to Gf, Tr[5*jT] and Tr[T] 
are real and the order of T is larger than 2. As discussed in the previous section, although 
not mandatory, they are general and reasonable assumptions on residual symmetries. It is 
interesting to see that phenomenologically viable mixing patterns they can lead to include 
only several well-known ones. This result provides a different way to understand those 
mixing patterns as natural consequences of discrete flavor symmetries. On the other hand, 
the no-go result obtained above implies that, to accommodate the mixing data including 
sin #i3, which is small but far from vanishing, some of the assumptions adopted in this section 
may have to be relaxed, which will be discussed in the next section. 



IV. BEYOND THE MINIMAL SCENARIO 



In this section, we focus on the consequences of allowing Tr [S^T] or Tr[Wj to be complex 
valued. In this case, the eigenvalues of Wi satisfying \ Pi = 1 can be written as e 2nmi ^ Pi , 
e 2mrn '/p* a nd e - 27 ™K+ m *)/p* . Without loss of generality, we require that < rii, rrii < Pi, 
rii + rrii 7^ and gcd(rii,pi) = 1. The solutions to Eq. (TT0|) and (ITT]) are given by 



(U 



PMNS) 



sm 



2 kir_\ 1 



m ) 



cos 



— 2tt 

2pi 



cos 



m* rii + mi 

27T COS 

2 P ,, 



2tt 



PMNS) jH 



_ . 2 kn kir 
2 sm — cos — 
m m 



-1 



cos 



2m 2p, t 



cos 



PMNS, 



k m; 
— + — 
2m 2pi 

krc 



2tt 
kir 



cos 
-1 



k 

2m 



rii + rrii 



cos 



2 sin" — cos — 

m m 



I 2tt 

2m 2p, 



cos 



cos 



2p 

k rii 
2m 2pi 
k rii + rrii 
2m 2pi 



2tt 



2tt 



2tt 



(20) 



(21) 



(22) 



11 



One can verify that Eq. (1151) and (fl~6]) are recovered when n« or m.; vanishes or n« + rrii = Pi. 

If all Si G Gf, as in the previous section, the unitarity condition must be respected and 
then from Eq. one has 

3 

^2 v i + 2 cos "fcm = -1, Vi = cos a„ iPi + cos a miPi + cos {a UtPi + a TO<p< ) (23) 



i=i 



and 



i=i 



where ctkm = 2irk/m, a niPi = 2irni/pi, etc. Unlike the minimal scenario where Tr[Wj] are 
real, the above two equations are much more involved and the complete set of solutions can 
hardly be obtained in this case. But, as in the minimal scenario, solutions for reasonably 
large pi can be exhausted by numerical calculations. Besides that, they are also useful for 
one to approach possible solutions quickly, especially when some Pi are given. 

As an example, let m = 4 and p\ = 12. Since cos (27mi/12) = ±V3/2, V x and V[ are 
expected to be a quadratic number 6 of the form a + by/3 except for particular choices of n\ 
and mi (e.g. n\ = 1 and rti\ = 5 or 6 which lead to | (Upmns)h I = ^)- Therefore, to satisfy 
Eq.( )23|) . it is reasonable to require that P2 or p 3 also gives rise to some quadratic terms with 
the factor \f?>. The simplest choice is to let, e.g. ps = 12. Then, without much effort, we 
find that setting p2 = 3 can lead to the mixing matrix 

/ ^3 + 73 y/2 

WUpmnsW = V2 2 y/2 

which presumably can be realized by the group A (432). 

Although the resulting sin 2 9 13 = 0.158 does not agree with the data well, the point shown 
by this example is that if any pi grows large, especially when the value of cos(27rnj/pi) or 
sin(27rnj/pi) is not a rational number 7 , then for generic n« and m«, it often occurs that 



6 By quadratic number we mean a number satisfying a quadratic equations with rational coefficients. 

7 A more precise condition can be given in terms of the algebraic degree of cos(2irrii / pi) or sin(27rrii/pi), 
which is the lowest degree of the algebraic equations with rational coefficients it satisfies. Except for 
m = 1 or 2, the algebraic degree of cos(27rn/m) is given by ip(m)/2 61] where (p(m) is Euler's ^-function 



or Euler's totient function 



62[. For example, for m = 3, 4, 6, ip(m)/2 — 1 and hence cos(27m/m) is a 
rational number of degree one. For m = 5, 8, 10 and 12, <p(m)/2 = 2 and cos(27rn/r7i) is a quadratic 
number of degree two. The degrees of numbers involved in algebraic operations, including addition, 
product, division, etc. can provide valuable information about the outcomes of these operations. 
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another pi have to acquire the same value. This also explains to some extend why identical 
columns up to permutations often show up in full-mixing matrices, as the one given above. 
Moreover, if m takes a value leading to a non-rational cos(27r/c/m), e.g. m = 5, 7, etc. the 
above unitary conditions are even harder to be satisfied since m appears only in one of them, 
i.e. eq. (1231) . This implies that even a large group can lead to only a rather limited number 
of mixing patterns. A slight improvement can be achieved by allowing Tr[T] to be complex, 
which will be discussed later in this section. 

As another interesting application, consider m = 3 corresponding to G e = Z3 which 
occurs frequently. When m = 3 one finds that Eqs. (1201 - 122"]) can be written as 



PMNS, 



,2 4 
= - COS 




cos 




COS 


1 3 


. 2 P P i . 




. 2 Ppi 





2p 



i" 



(25) 



where p = a, (3 or 7, p 



Pi, 



n. 



ft; 



mj and 
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in which and p 7 j are chosen to be the smallest (positive) integers satisfying the above 
equations and H(x) is the unit step function defined as H(x) = 1 (for x > 0) and H(x) = 
(for x < 0). As shown by Eq.( ]2"5"]) . the similarity among the expressions for | {Upmns)™ 1 2 
indicates that and p 7i can also play the role of Pj. Therefore, by assigning pj, p^j and 
p 7 j to the orders of S{T one can construct a mixing matrix with identical columns and rows 
up to permutations. For example, {p 3 , n 3 , m 3 , m, k} = {5, 2, 3, 3, 1} results in p^ = p 7i = 15 
and the squared mixing vector (6 — 2\/5, 3 + \/5, 3 + \/5j / 12. Then one may let pi_2 = 15 
and the mixing matrix 

v^+l v^+l VW-V2~\ 

i 

v/TO-v 7 ^ v^+1 v^+l 
v v^+1 V10-V2 V5 + 1 j 

can be obtained by permuting the elements in the third column. A possible issue of the 
mixing patterns alike is that to fit the experimental data somewhat large corrections are 
required since in every column there is an element equal to | (Upmns)^ I- Also note that in 

13 



\U 



PMNS\ 



2V6 



this example TrpHy is real but TrfWi^] are not. In general, to produce mixing matrices like 
the one above, some Tr[Wj] must be complex valued, otherwise the matrix ||{7pmats|| should 
have two identical rows. 

We now turn to phenomenological implications. As we know, most known full-mixing 
patterns cannot fit the mixing data exactly within the experimentally allowed range. Nev- 
ertheless, by using the method discussed above, it is not hard to find one. For example, 
{PiiP2,P3, m } = {20, 3, 10, 3} can lead to: 

/ \/9 + a/5 + c 2a/2 a/7 - a/5- c\ 



\u\ 



2a/6 



6-2a/5 2y/2 V10 + 2V5 



, a/9 + a/5 - C 2a/2 a/7- a/5 + C 



where c = y 6(5 — a/5). The mixing angles can be extracted: sin 2 # 13 = 0.029, sin 2 6*23 = 0.62 
and sin 2 12 = 0.34. The obvious problem with this example is that it requires a large group, 
presumably A(600). In fact, this result is confirmed by group theoretical calculation in 
a recent work 54j] in which a scan of groups with orders up to 1536 is performed. Our 
calculation does not need any prior knowledge about A(600) and, as discussed above, it is 
the solution that sug gests A(600) as a candidate. Besides A(600), the other two mixing 
patterns found in 54| for A(1536) and (Z 18 x Z 6 ) x S3 can also be obtained from Eqs. (I20H22I) . 
Just in case it is needed, below we give the exact value of the mixing matrix for A (1536) 

( \j2 + a/2 + a/I a/2 a/2 - a/2 + a/I ^\ 
v^^T a/2 v/2T^I 
a/2^7 V2 



11711 = 4* 
a/6 



V 



v / 2T^ J 



where c± = y 2 ± y 2 =p a/3. The numerical values of sin 2 9{j can be found in [54]. As to 
the group (Zi§ x Z$) x S3, the corresponding mixing matrix does not admit an expression 
in terms of exact values as the one shown above 8 . 



Using algebraic number theory, one can show that mixing elements can be expressed in terms of exact 
values only when some conditions are satisfied. For example, in the case where m = 3, the orders of SiT 
must be products of Fermat numbers. By exact values we mean numbers that can be expressed in terms 
of rational numbers involving only sums, products and square roots, such as (2 — a/3) /3. This also implies 
that among all the mixing patterns that can be expressed in exact forms, only very particular ones such 
as those discussed in this paper, can be produced by discrete symmetries. Therefore, one can show that 
some mixing patterns, e.g. the hexagonal mixing, cannot be produced by discrete symmetries at least in 
the framework discussed in this work. More detailed discussion will be given elsewhere. 



14 



Besides [541 ] . it is also reported in [53| that no SU(3) subgroup of order less than 512 
can produce the full mixing data. For us this is almost an expected result since from 
Eqs. (I20H22P it is easy to see that small sin# 13 requires one or more angles involved to be 
close to it/ 2 or 3it/2 and hence, together with the reality conditions, it would require Pi 
or m to be large. Eqs. ( I201 - r2"2"j) can also be evaluated numerically and the result agrees 
with them, as expected. Since in general mixing parameters also receive contributions from 
higher-order corrections including radiative corrections, this result implies that it might be 
more plausible to introduce sizable corrections (or free parameters) to the mixing patterns 
obtained from discrete symmetries, which can be considered as lea ding order contributions 



50 



52]. 



[431 ] . Partial-mixing is also a reasonable option, as discussed in 

Finally we comment on the case where Im(Tr[T]) ^ 0. When it occurs, the unitarity 
conditions given by Eq. (l2"3"|) and (124"]) . which follow from fl9]), should be modified to 



^2 v i = 0) Vi = cos a n iPi + cos a mm + cos («„,,,, + a mm ) (26) 



i=i 

and 



3 

X " V- = V' A , V- = sin a niPt + sin a m%n - sin (a niP , + a mm ) . (27) 



Although more complicated and even harder to be satisfied, it allows new mixing patterns 
for particular pi. As an example, let m = p^ = 7, n 4 = 1. The quickest way to find a solution 
in this case is to use the trick of assigning the value of p^ or m to another pi, which works 
again fortunately. After that, the remaining pi can be easily derived. One can verify that a 



solution can be found and it leads to another PSL(2, Z-j) mixing matrix given in 45] (see 
Eq. (40) there). 

Although a solution is found in the example above, it is not hard to see that it is very 
likely that this solution is the only solution for the case where m = 7, which shows that if m 
is associated with a cos(27r/m) with a high degree (see footnote 7), the unitary conditions 
are hard to be satisfied and solutions may not always exist. It also provides a further support 
to our finding that using large groups may not help too much for the purpose of producing 
full-mixing matrices that can fit all the experimental data. 
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V. CONCLUSIONS 



As shown in the previous sections, the method discussed in this paper can provide an 
efficient way to derive neutrino mixing patterns determined by underlying discrete flavor 
symmetries. Although most of our attention is paid to full-mixing patterns, the formalism 
developed in this paper can also be used to calculate mixing elements in the partial-mixing 
case. It can be used to cross check with group theoretical calculations or as a simpler way 
to reach some known results. It also gives the necessary conditions that can lead to the 
flavor symmetry corresponding to a mixing matrix or mixing-matrix column derived by this 
method. Since the properties of most SU (3) subgroups can be found in existing literatures, 



e.g 



63 



-|65|. one may find an appropriate group without much trouble in most cases and 
the physical viability can be verified concretely, as shown by the A (600) example. 

Moreover, it can also shed new light on the relation between neutrino mixing and discrete 
symmetries. We find that in the minimal scenario where some general assumptions are 
adopted, except phenomenologically disfavored mixing patterns, discrete symmetries cannot 
lead to mixing patterns other than the well-known TBM, BM and golden ratio mixings. 
This result not only provides a new way to understand these mixing patterns as simple and 
natural consequences of discrete symmetries, but also indicates clearly the limitation of the 
minimal scenario and how to go beyond it. Another interesting result can be found by this 
method is that, although relaxing some assumptions can lead to new mixing patterns, to 
fit the mixing data, relative large groups are required for full-mixing matrices unless sizable 
corrections are allowed. 
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